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Abstract 

We implement, at next-to-next-to-leading order, the QCD corrections to Standard 
Model Higgs boson production in association with vector bosons at hadron colliders, 
qq HV with V = W, Z . They consist of the two-loop corrections to the Drell-Yan 
process for the production of off-shell vector bosons, qq — > V*, and in the case of Z 
final states, of the additional contribution from heavy-quark loop mediated processes, 
in particular gg — > HZ. For the Higgs boson masses relevant at the Tevatron and 
the LHC, ^ 200-300 GeV, the two-loop corrections are small, increasing the pro- 
duction cross sections by less than 5% and 10%, respectively; the scale dependence is 
reduced to a level of less than a few per cent. This places these processes among the 
most theoretically clean Higgs boson production channels at hadron colliders. 
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1. Introduction 



The Standard Model (SM) predicts the existence of a scalar particle, the Higgs boson, that 
is the remnant of the electroweak symmetry-breaking mechanism that generates the weak 
gauge boson and fermion masses |Tj. The search for this particle is the primary mission of 
present and future high-energy colliders. The Higgs boson can be discovered at Run II of the 
Tevatron if it is relatively light, < 200 GeV, as suggested by the fits to the high-precision 
electroweak data j2], and if sufficient integrated luminosity is collected [21111 • Higgs bosons 
with masses up to ~ 1 TeV, a value beyond which perturbation theory is jeopardized 
in the SM, can be probed at the upcoming LHC jHllll. 

One of the most important Higgs boson production mechanisms at hadron colliders is 
the Higgs-strahlung process, i.e. the associated production of Higgs and weak gauge bosons, 
qq HV with V = W, Z At the Tevatron, Higgs particles can be mainly produced 
in the channel qq — HW with the W boson decaying into iu pairs [with i = e, fi] and the 
Higgs boson decaying into bb or W^W~ pairs [71. At the LHC, a plethora of production 
channels can be used to search for the Higgs particle; one of the principal detection modes is 
expected to be the gluon-gluon fusion process, gg ^ H [HUHl, with the signatures i7 — > 77 
OT H ^ ZZ^*\ WW M in, respectively, the low and high Higgs mass ranges. However, 
although with more difficulty, the Higgs boson can also be detected through the channels 
qq — > HW/HZ, in particular in the 77 plus lepton final states [TUllII]. These processes 
could play a very important role in the determination of the Higgs boson properties |12j . 

It is well known that, in hadronic collisions, the lowest-order (LO) cross sections are 
affected by large uncertainties arising from higher-order QCD corrections. If at least the 
next-to-leading order (NLO) radiative corrections are included, the cross sections can be 
defined properly and their unphysical variation with the scales are stabilized. To have an 
even better control on the theoretical prediction, the next-to-next-to-leading order (NNLO) 
corrections, which are in general very complicated to calculate, are desirable. Up to now, 
the NNLO corrections to SM Higgs boson production at hadron colliders are known only for 
the gg ^ H mechanism [^ in the infinite top quark mass limit ^. 

In this paper, we will discuss the NNLO order, i.e. the 0{a1), corrections to the pp — * HV 
production cross sections [hereafter, we will use the notation pp for both pp and pp\ . Part 
of these corrections are simply those of the Drell-Yan process pp — > V\ however, in the case 
of pp HZ, additional corrections are due to diagrams involving the Zgg vertex as well as 
to loop-induced gg fusion diagrams. We show that while the NNLO corrections increase the 
cross sections by only 5 to 10%, the scale dependence of the latter is drastically reduced, 
making these channels among the theoretically cleanest Higgs production mechanisms. 

^The NNLO corrections to the bottom-quark fusion mechanism, bb H + X, which plays an important 
role in supersymmetric extensions of the SM, have also been derived recently 
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The paper is organized as follows. In the next section, we describe the known behaviour 
of the production cross section at leading and next-to-leading orders in QCD. In section 3, we 
summarize the main contributions to the cross section at NNLO: the Drell-Yan corrections 
and the additional corrections to HZ production due to gluon-gluon-initiated processes. The 
numerical results for the -R'-factors, the variation with the renormalization and factorization 
scales are summarized in section 4. A short conclusion and some remarks on the remaining 
uncertainties on the production cross section are given in the last section. 

2. LO and NLO cross sections 

The associated production of Higgs and gauge bosons, Fig. 1, is one of the simplest produc- 
tion mechanisms at hadron colliders: the final state does not feel strong interactions, which 
affect only the quark and antiquark initial state. In fact, this process can be viewed simply 
as the Drell-Yan production of a virtual W or Z boson, which then splits into a real vector 
boson and a Higgs particle. Denoting by k the momentum of the virtual gauge boson, the 
energy distribution of the full subprocess can be written at leading order as 

^{qq - HV) = a{qq - V*) x ^(V^* - HV) , (1) 

where, in terms oi < < = s with s the centre-of-mass energy of the subprocess and 
the usual two-body phase-space function \{x,y; z) =(1 — xjz — yjzY — ^xyjz^^ one has 

^iV^ HV) = GFM'y\'l\Ml,Ml-k^) ( X{Ml,Ml-k^) \ 

dP^ ^ 2V27r2 (F-M2)2 12M2/F J' ^' 




V = w.z 



Q H 
Figure 1: Diagram for associated Higgs and vector boson production in hadronic collisions. 

The total cross section for the subprocess is obtained by integrating over k"^ : 



^Lo(gg - VH) = + al)Xy\M^y, Ml- s) 

ZooTTS 



A(M^,M^;g) + 12M^/g 
(1 - Ml/sf 



(3) 



where the reduced quark couplings to gauge bosons are given in terms of the electric charge 
and the weak isospin of the fermion as: = 2/^ , Vq = 21^ — 4:Qqs'^ ior V = Z and 



V2ioTV = W, with sL = l-c 



sin^ 9iv. 
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The total hadronic cross section is then obtained by convoluting eq. (3) with the parton 
densities and summing over the contributing partons 



cThoipp VH) 



dr > — — (Tlo(s 



T"0 



dr 



TS 



(4) 



where Tq = {My + Mfj^/s = Mfjy/s, s being the total hadronic cm. energy, and the 
luminosity is defined in terms of the parton densities defined at a factorization scale fip- 

In fact, the factorization of the pp HV cross section in eq. (1) holds in principle at 
any order of perturbation theory in the strong interaction and we can thus write^ 



da 



{pp HV + X) = a{pp 



V* + X)x —{V*^HV) 



(5) 



where dF/d/c^ is given by eq. (2). Therefore, the QCD corrections to the Higgs-strahlung 
process, derived at NLO in Refs. P^lSj . are simply the corrections to the Drell-Yan process 
[inillZ], as pointed out in Ref. [TT] . 

At NLO, the QCD corrections to the Drell-Yan process consist of virtual corrections 
with gluon exchange in the qq vertex and quark self-energy corrections, which have to be 
multiplied by the tree-level term, and the emission of an additional gluon, the sum of which 
has to be squared and added to the corrected tree- level term; see Fig. 2. 
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1/ 1/ 1. 

Figure 2: NLO QCD corrections to the vector boson-quark-antiquark vertex. 

Including these contributions, and taking into account the virtuality of the vector boson, 
the LO cross section is modified in the following way 



with 



Aa. 



QQ 



An 



Q9 



c"nlo 



TT 
IT 



0"LO 



A 



QQ 



An 



dr 



to/t 



QQ ' 



dz (Jlo(t2;s) UJgg{z) , 



(6) 



drj: 



Q,Q 



dO^ 
dr 



t"o/t 



dz a^oiTzs) ujgg{z) 



(7) 



^This is only valid at first order in the electroweak coupling; at two-loop order in G_f, QCD corrections 
to the final state should also be taken into account. In addition, for the process pp — > HZ, an additional 
contribution will appear at 0{a'^) as will be discussed later. 



4 



with the coefficient functions 

,,2 



vr 



__4)5(l_^) + 2(l + ^2^ 



log(l - z) 
l-z 



--Pggiz) log 



1 — z)'^TS 



1 

+ - 



1 + 62 - 7z' 



where nn denotes the renormahzation scale and Pqq, Pgg are the well-known Altarelli-Parisi 
splitting functions, which are given by 



1 + z^ 

L(i-^)+ 



-m~z) 



1 r 
2 



z' + {l- zf 



(9) 



The index + denotes the usual distribution F+{z) = F{z) - 5(1 - z) J^dz'F{z'). Note 
that the cross section depends explicitly on log(/i|,/(5^); the factorization scale choice = 
= Mfjy therefore avoids the occurrence of these potentially large logarithms. The renor- 
mahzation scale dependence enters in the argument of and is rather weak. In most of 
our discussion, we will set the two scales at fip = fiji = Mhv- For this choice, the NLO 
corrections increase the LO cross section by approximately 30%. 



3. The NNLO corrections 



V* 



The NNLO corrections, i.e. the contributions at 0{a'j.), to the Drell-Yan process pp 
consist of the following set of radiative corrections [see also Fig. 3a-c]: 

a) Two-loop corrections to qq ^ V*, which have to be multiplied by the Born term. 

b) One-loop corrections to the processes qg qV* and qq gV*, which have to be 
multiplied by the tree- level gq and qq terms initiated by the diagrams shown in Fig. 2. 

c) Tree-level contributions from qq,qq,qg, gg V*+ 2 partons in all possible ways; the 
sums of these diagrams for a given initial and final state have to be squared and added. 





b) c) 

Figure 3: Diagrams for the NNLO QCD corrections to the process qq — ^ W* . 
These corrections have been calculated a decade ago in Ref . [Ej and recently updated (TH] . 
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However, these calculations are not sufficient to obtain a full NNLO prediction: the cross 
sections must be folded with the NNLO-evolved parton distribution funtions (PDFs), which 
are necessary. The latter require the calculation of Altarelli-Parisi splitting functions up to 
three loops and, to this day, the latter are not completely known at this order. Nevertheless, 
a large number of moments of these functions are available [201; when these are combined 
with additional information on the behaviour at small x, we can obtain an approximation 
of the splitting functions at the required order. The NNLO MRST |25 parton distributions 
follow this approach and can therefore be adopted for this calculation. 

In the case of pp HZ production, because the final state is electrically neutral, two 
additional sets of corrections need to be considered at O^a"^). 

Contrary to charged W bosons, the neutral Z bosons can be produced via an effective 
Z-gluon-gluon coupling induced by quark loops (Fig. 4). This can occur at the two-loop 
level in a box+triangle diagram in qq Z* [which has to be multiplied by the Born term], 
or at the one-loop level where vertex diagrams appear for the qq — > gZ* and qg — qZ* 
processes [and which have to be multiplied by the respective 0{as) tree-level terms]. 




Figure 4: Diagrams for the QCD corrections to qq — > Z* not present in qq W* . 

Because gluons have only vector couplings to quarks and the effective Zgg coupling must 
be a colour singlet, only the axial-vector part of the Zqq coupling will contribute as a 
consequence of Furry's theorem. The axial-vector of the Z coupling to quarks, Og = 2/^, 
differs only by a sign for isospin up-type and down-type quarks, so that their contribution 
should vanish in the case of quarks that are degenerate in mass. Thus, in the SM, only the 
top and bottom quarks will contribute to these topologies^. 

These corrections have been evaluated in Refs. P^ l^ and have been shown to be ex- 
tremely small. The one-loop corrections give relative contributions that are less than a few 
times 10~^ and thus completely negligible. The two-loop contribution is somewhat larger. 
However, at Tevatron energies and for partonic cm. energies close to the Z boson mass, 
the contribution is still below the 1% level, and it is even smaller at the LHC [221 • These 
corrections can therefore be safely neglected. 

•^Note that additional contributions involving heavy top quarks in gluonic two-point functions are also 
present; however, they give zero contribution when the top quark is decoupled. 
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Another set of diagrams that contribute at 0{a'^) to ZH and not to WH production 
[again because of charge conservation] is the gluon-gluon-initiated mechanism gg HZ 
P^I25j . It is mediated by quark loops [see Fig. 5] which enter in two ways. There is first a 
triangular diagram with gg ^ Z* ^ HZ, in which only the top and bottom quark contribu- 
tions are present since, the Z boson couples only axially to the internal quarks, because of 
C-invariance, the contribution of a mass-degenerate quark weak-iso doublet vanishes. There 
are also box diagrams where both the H and Z bosons are emitted from the internal quark 
lines and where only the contribution involving heavy quarks which couple strongly to the 
Higgs boson [the top quark and, to a lesser extent, the bottom quark] are important. 




Figure 5: Diagrams for the gg HZ process, which contributes to 0{a 



At the LHC, the contribution of this gluon-gluon fusion mechanism to the pp — > HZ total 
production cross section can be substantial. This is due to the fact that the suppression of 
the cross section by a power (as/vr)^ is partly compensated by the increased gluon luminosity 
at high-energies. In addition, the tree-level cross section for qq HZ drops for increasing 
cm. energy and/or Mh values, since it is mediated by s-channel gauge-boson exchange. 

We have recalculated the cross section for the process gg HZ at the LHC energy 
^/s = 14 TeV, taking into account the full top and bottom quark mass [rrib = 5 GeV and 
rrit = 175 GeV] dependence. The two contributing triangle and box amplitudes interfere 
destructively, as found in Ref. |2S]- Our results agree with those given in the figures of this 
article, once we take the same kinematical configuration, inputs and PDFs [^/s = 17 TeV 
and rrit = 80, 140 or 200 GeV]. The cross section for this process is of course negligible at 
the Tevatron, because of the low gluon luminosity and the reduced phase space. 

4. Numerical analysis 

The impact of higher order (HO) QCD corrections is usually quantified by calculating the 
/T-factor, which is defined as the ratio of the cross sections for the process at HO (NLO or 
NNLO), with the value of and the PDFs evaluated also at HO, over the cross section at 
LO, with as and the PDFs consistently evaluated also at LO: 

auo{pp^HV + X) 

^HO = 7 7777\ • (lOj 

o"Lo(PP HV) 
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A kind of /'^-factor for the LO cross section, -R'lo; can also be defined by evaluating the latter 
at given factorization and renormalization scales and normalizing to the LO cross sections 
evaluated at the central scale, which, in our case, is given hj fip = fiR = Mhv- 

The ii'-factors at NLO and NNLO are shown in Figs. 6 and 7 (solid black lines) for, 
repectively, the LHC and the Tevatron as a function of the Higgs mass for the process 
pp HW; they are practically the same^ for the process pp HZ when the contribution 
of the gg — > HZ component is not included. The scales have been fixed to fip = fJ^R = Mhv 
and the MRST sets of PDFs for each perturbative order are used in a consistent manner. 

The NLO i^-factor is practically constant at the LHC, increasing only from -R'nlo = 1-27 
for Mh = 110 GeV to Knlo = 1-29 for Mh = 300 GeV. The NNLO contributions increase 
the /T-factor by a mere 1% for the low Mh value and by 3.5% for the high value. At the 
Tevatron, the NLO /^-factor is somewhat higher than at the LHC, enhancing the cross 
section between i^NLO = 1.35 for Mh = 110 GeV and K^lo = 1.3 for Mh = 300 GeV with a 
monotonic decrease. The NNLO corrections increase the /T-factor uniformly by about 10%. 
Thus, these NNLO corrections are more important at the Tevatron than at the LHC. 

The bands around the i^-factors represent the variation of the cross sections when they 
are evaluated at renormalization and factorization scale values that are independently varied 
from ^Mhv ^ f^F (/^i?) < 3Mj^y, while the other is fixed to fiR (fip) = Mhv] the normal- 
ization provided by the production cross section evaluated at scales fJ^F = fJ'R = Mhv- As 
can be seen, except from the accidental cancellation of the scale dependence of the LO cross 
section at the LHC, the decrease of the scale variation is strong when going from LO to 
NLO and then to NNLO. For Mh = 120 GeV, the uncertainty from the scale choice at the 
LHC drops from 10% at LO, to 5% at NLO, and to 2% at NNLO. At the Tevatron and for 
the same Higgs boson mass, the scale uncertainty drops from 20% at LO, to 7% at NLO, 
and to 3% at NNLO. If this variation of the cross section with the two scales is taken as an 
indication of the uncertainties due to the not yet calculated higher order corrections, one 
concludes that once the NNLO contributions are included in the prediction, the cross section 
for the pp — > HV process is known at the rather accurate level of 2 to 3%. 

Finally, we present in Fig. 7 the total production cross sections at NNLO for the processes 
qq HW and HZ at both the Tevatron and the LHC as a function of Mh- In the case 
of the HZ process, the contribution of the gg ZH subprocess to the total cross section 
is not yet included, but it is displayed separately in the LHC case. For Higgs masses in the 
range 100 GeV < Mh ^ 250 GeV, where a{qq — > HZ) is significant, cr(gg HZ) is at the 
level of 0.1 to 0.01 pb and represents about 10% of the total cross section for low Mh- 

^Because of the slightly difTcrcnt phase space and scale, the if -factor for pp ZH is not identical to the 
-ftT-factor for pp WH . However, since (Mf — M^)/s is small and the dependence of dF in eq. (3) on 
is not very strong in the range that we are considering, the if-factors for the two processes are very similar. 
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NLO 



^ 1.4 
U 

ffil.35 

-X 1.3 
^1.25 
1.2 
1.15 
1.1 
1.05 
1 

0.95 

120 140 160 180 200 220 240 260 280 300 

Figure 6: The K-factors for pp — * HW at the LHC as a function of Mh at LO, NLO 
and NNLO (solid black lines). The hands represent the spread of the cross section when the 
renormalization and factorization scales are varied in the range \Mhv < yLtfl {iJ'f) < SMj^y, 
the other scale being fixed at fipif^R) = Mhv- 
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Figure 7: The same as in Fig. 6 hut for the Tevatron case, pp HW . 



9 



Its relative magnitude increases for higher Higgs masses; for very large Mh values, it 
reaches the level of the Drell-Yan cross section. However, for these large Higgs masses, the 
total production rate is too small to be useful in practice. Note that at the Tevatron, as 
expected, the cross section of the gg HZ subprocess is very small, barely reaching the 
level of (y[gg HZ) ~ 0.2 fb for Mh = 120 GeV. The contribution of this subprocess can 
be safely neglected in this case. 




100 150 200 250 300 100 150 200 250 300 



Mh [GeV] Mh [GeV] 

Figure 8: The total production cross sections at NNLO forpp HW and HZ at the LHC 
(left) and the Tevatron (right) as a function of Mh- The MRST parton densities have been 
used. The contribution of the gg HZ process is shown separately in the case of the LHC. 

5. Concluding remarks 

We have discussed the Higgs-strahlung processes pp HV with V = W, Z ai NNLO in 
strong interactions. We have shown that the Drell-Yan-type corrections increase the total 
production cross sections by up to 3% at the LHC and by up to 10% at the Tevatron in 
the Higgs mass range relevant at these colliders. Because of the larger gluon luminosity at 
high energies, the additional contribution due to the gg HZ subprocess can be relatively 
important at the LHC, and for Higgs boson masses far above 200 GeV, it becomes comparable 
with the qq — * HZ production cross section, but the total rate is then rather small. The 
scale dependence is strongly reduced from NLO, where it can reach a level close to 10%, to 
less than 2-3% at NNLO. These NNLO corrections are of the same order, but of opposite 
sign, as the 0{a) electroweak corrections to these processes, which have been calculated very 
recently [21]. 
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Together with the effects of higher-order corrections, the uncertainties due to the PDFs 
dominate the theoretical error on the production cross section. Recently, the CTEQ j2ZI and 
MRST j2H| collaborations introduced a new scheme, which gives the possibility of controlling 
these errors: in addition to the nominal best fit PDFs, they provided a set of 2N PDFs 
[N = 20 for CTEQ and 15 for MRST] at NLO [they are not yet available at NNLO], 
corresponding to the minima and maxima of the eigenvectors of the matrix error of 
the fitting parameters. Adding the maximum and minimum deviation for each eigenvector 
in quadrature, we obtain an error on the total cross section at NLO of less than 5% for 
Mh < 300 (200) GeV for LHC (Tevatron) energies 

All these features make Higgs-strahlung one of the most theoretically clean Higgs boson 
production channels at hadron colliders. This will be of great importance when it comes 
to the determination of the properties of the Higgs boson and the measurement of its cou- 
plings. The systematic uncertainties originating from higher-order corrections and structure 
functions being small [a better determination of the parton distribution functions can be per- 
formed in the future], the Higgs-strahlung process will provide a clean determination of the 
HVV couplings times the Higgs branching ratios [the latter being measured in other Higgs 
production processes] if enough integrated luminosity is collected to make the statistical 
errors also small^. 

This analysis can straightforwardly be extended to the case of the Minimal Supersym- 
metric extension of the SM p. The two CP-even Higgs bosons, $ = /i and H, can be 
produced in the same channels, pp and the LO cross sections are the same as the 

one for the SM Higgs boson, except that they are suppressed by global coupling factors, 
< 

9^vv — ^- '^^^ standard QCD corrections are again similar to those discussed here; 
the additional corrections due to supersymmetric partners of quarks and gluons have to be 
added [note that the triangular Zgg diagrams with scalar loops give zero contribution here]. 
This analysis can even be extended to the case of the associated production of CP-even and 
CP-odd Higgs particles, pp ^A. Here again, the cross section can be factorized into the 
product of pp — > V* production and V* decays, and the QCD corrections are thus the 

same as those discussed here. A small additional contribution originating from the one-loop 
subprocess gg ^A ^Tj has to be added, though. 

Acknowledgements: We thank Samir Ferrag and Michael Spira for discussions. 

^Note that an additional systematic error of about 5% arises from the pp luminosity. To reduce all these 
uncertainties at the LHC, it has been suggested (SUj to use the Drell-Yan processes pp ^ W,Z with the 
subsequent leptonic decays of the gauge bosons as a means of measuring directly the quark and antiquark 
luminosities at hadron colliders. The errors on the cross sections for all processes dominated by qq scattering, 
such as the Higgs-strahlung discussed here, when normalized to the Drell-Yan rate, would lead to a total 
systematical uncertainty of less than 1%. In this case, the dominant part of the if-factor for Higgs-strahlung 
will drop out in the ratio. 
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